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A decoupled system of hyperbolic partial differential equations for linear perturbations around any 
spatially flat FRW universe is obtained for a wide class of perturbations. The considered perturbing 
energy momentum-tensors can be expressed as the sum of the perturbation of a minimally coupled 
scalar field plus an arbitrary (weak) energy-momentum tensor which is covariantly conserved with 
respect to the background. The key ingredient in obtaining the decoupling of the equations is 
the introduction of a new covariant gauge which plays a similar role as harmonic gauge does for 
perturbations around Minkowski space-time. The case of universes satysfying a linear equation of 
state is discussed in particular, and closed analytic expressions for the retarded Green's functions 
solving the de Sitter, dust and radiation dominated cases are given. 
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1. Introduction 

The theory of linearized perturbations around cosmo- 
logical backgrounds is a cornerstone in the understanding 
of structure formation in the universe, because it provides 
the connection between what we see today, and the ini- 
tial conditions in the early universe Q. This theory has 
a long history, beginning (in the context of general rel- 
ativity) with the pioneering work of E. Lifshitz gj. An 
important step was the introduction by J. Bardeen of 
gauge invariant potentials to describe the metric pertur- 
bations II . This is the so called standard treatment of 
cosmological perturbations, and many authors have elab- 
orated this approach and used it to deal with the analysis 
of observations [|| . The standard treatment is based on 
the decomposition of the metric perturbations in scalar, 
vector and tensor parts according to their behaviour un- 
der the rotation group of isometries of the FRW back- 
ground fl. In the standard approach these three kinds 
of perturbations are treated independently, and most at- 
tention has been devoted to the discussion of scalar per- 
turbations in the particular case in which the two gauge 
invariant scalar potentials coincide. This particular case 
is well motivated because it arises naturally when dealing 
with hydrodynamic or scalar field perturbations. In these 
cases, the space-space components of the perturbing en- 
ergy momentum-tensor take the form ST? oc S 8?, where 
S is a scalar, and in turn this implies the equality $ = '!', 
between the two gauge invariant Bardeen potentials for 
the scalar sector. An important technical consequence 
of this equality is that one has to deal only with a sec- 
ond order differential equation for <3>; in contrast with the 
general case of scalar perturbations, in which the stan- 



dard formalism yields two coupled equations for $ and 

In this letter we generalize previous results that we first 
obtained for a dust filled universe J(| , and present a new 
approach to cosmological perturbations, which allows a 
unified treatment of all types of perturbations: scalar, 
vector and tensor. The method follows a very similar 
path to that used for computing the gravitational ra- 
diation emitted by astrophysical sources in Minkowski 
space-time, using harmonic gauge As a result of our 
approach, we obtain a decoupled system of hyperbolic 
partial differential equations for all components of the 
metric perturbation /i^„. This decoupling of equations 
holds for a very wide class of perturbations, which we 
define below, including the presence of seeds like cos- 
mic strings or primordial black holes |?J . In addition the 
equations for perturbations are formulated in a covari- 
ant way, and in our opinion they are best fitted to study 
linearized quantum gravity on spatially flat cosmological 
backgrounds. 



2. The background model 

We choose a spatially flat FRW universe as back- 
ground. This simplifies matters a little, at not a very 
high cost concerning physics, because the recent analysis 
of observations of supernova type la || and of acous- 
tic peaks in the distribution of CMB temperature versus 
angular momentum JIo[] , favour a flat universe. In addi- 
tion, inflationary models of the very early universe also 
favour flat universes |pd| . The metric of such background 
model written in conformal time r\ and spatial cartesian 
coordinates x, reads W3 
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ds 2 = 7m „ dx^dx" = a(r]) 2 (-drj 2 + dx 2 ) , (1) 
where the scale factor a(rj) is an arbitrary function of 



2 



time. Thus, the equations for perturbations that we are 
going to derive will apply to any flat FRW background, 
and in particular to multifluid universes, like our own 
universe. To describe the matter producing this back- 
ground, we shall use a scalar field <f> minimally coupled 
to gravity, i.e. we shall consider the action 



S [gap, 



lkI d4xV ~ gR 

d A x^-g{^g a ^ ia fa +V(4>)) (2) 



The full metric g a p and the scalar field <fi are to be 
developed up to first order in perturbations, in the form 
9a.p = la/3 + h a p, and (f> = 4> + 6<f>. The zero order parts 
describe the background. It is important to emphasize 
that this description of the background, by means of the 
scalar field 0, does not place any restriction on the scale 
factor a(rj). Any desired function a(r]) can be obtained 
by an appropriate choice of the potential V((f>). In this 
sense, the scalar field 4> is not necessarily to be thought 
of as a fundamental field filling the universe, but as a 
convenient way of parametrizing the matter content of 
the universe at cosmological scales; as for example cold 
matter plus a cosmological constant, or whatever model 
future observations will support. 

Using this scalar field parametrization of the matter, 
the zeroth order energy-momentum tensor can be written 
in the perfect fluid covariant form 



8" + U a U^(j), a <f>,p 



(3) 



where u v is the background velocity field of the fluid. 
The values of the background pressure and density can 
be derived from (^), and the covariant conservation of 
T„, with respect to the background metric, yields the 
zeroth order equation of motion for the field <f> 



I" = a' 1 <K 



(7) 



Writting a = exp Q, the Hubble parameter is given by 



H = — = tl exp(-fi) , 
the Einstein equations for the background are 



\2 



o + - n 2 = -47rg 



4 -a 2 V{$) 



(8) 



(9) 



(10) 



and the integrability condition (H) reads 



+ 2f2</> + a : 



dV(<£) 



Q 



(11) 



For example, in the particular case of a flat FRW back- 
ground with a linear equation of state p = a p, the cor- 
responding potential for the scalar field is 



v(4>) 



3(1- a) Hi 
16tt£ 



exp (-2^/6(1 + a)7rg (<j> - <j> )) 



(12) 



and the solutions of the Einstein equations for the Hub- 
ble parameter, the background scalar field (j>, and the 
conformal Hubble parameter O, are 



3a 



H( V )=H ^ 



Vo \ 1 + 3a 



(13) 



dV(4>) 



(4) 



where the vertical bar | denotes the covariant derivative 
with respect to the background metric. 

In coordinates (r),af), and using a dot for derivatives 
with respect to conformal time ry, the background values 
of pressure, density, and velocity take the form: 



P=\a- 2 4> 2 + V($) , 



(5) 
(6) 



and 



= 90 



1 / 3 + 3a rj_ 
l + 3aV 2tt£ U r]o 



and 



(l + 3a)r? 



(14) 



(15) 



the label refers, as usual, to the present epoch values 
of the corresponding quantities, and the evolution of the 
density follows the critical value law p = 3H 2 /8ttQ. 



3. Structure of perturbations and gauge invari- 
ance 

Now we are going to specify the class of perturbations 
that we will allow on the background. The equations 
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for perturbations come out, as usual, by developing the 
Einstein tensor up to first order in perturbations G^ = 
+ 8G", and splitting the Einstein equations into a 
zeroth order part 

G/ = -8^T/ , (16) 

which describes how the background is created by T^, 
and a first order part 

JG/ = -8tt0£T/ , (17) 

which gives the equation governing perturbations. 

Expansion of the Einstein tensor up to first order in 
pcturbations yields the covariant expression 



background metric. On the other hand, what we call 
the free perturbation corresponds to additional matter 
(more or less exotic) that can exist in the universe in ad- 
dition to the main component, and which moves as test 
matter on the background, i. e. it fulfills equation (pl|). 
This extra matter could be, for instance, topological de- 
fects as cosmic strings or domain walls produced by phase 
transitions in the early universe ||l4|| , or any other seed 
perturbation as for example primordial black holes. But 
also, ordinary baryonic matter, amounting only to two or 
three percent of the total energy density in the present 
universe, can be considered as a free perturbation too. In 
the last case the intrinsic perturbation 8T^ U would cor- 
respond entirely to nonbaryonic dark matter or energy. 

Now the expression for ST^'" is obtained from the 
action (En, and in turn, it can be decomposed as 



2<5G/ 



A \v 
/i I A 



R a 0lfj 



v |A 
a/3 



1 



Rib 



(18) 



where ip£ = h£ — ^8" h is the so called trace reversed 
graviton field. On the other hand the structure of the 
perturbing energy-momentum tensor ST^ v is constrained 
by the Bianchi identity. Expanding the full Bianchi iden- 
tity to first order in perturbations, and taking into ac- 
count the covariant conservation of TV" with respect to 
the background metric, one obtains the constraint 



6T/ W + 8T» a T^ 



ST a T v = 



(19) 



where 8T* is the first order part of the full metric con- 
nection. The inhomogeneous equation (|l^) must be ful- 
filled by any perturbing energy momentum tensor. Thus, 
any perturbing 8T^ can be decomposed in the form 



5T u = 6T (I), . 



s T (F)v 



(20) 



where 8Tu is a particular solution of the constraint 



equation (|19|), and ST^ 
neous equation 



(F)v 



ST, 



(21) 



In addition, we will assume that ST^ V does not func- 
tionally depend on the metric perturbation h^ v . We shall 

call STp" the intrinsic perturbation, and 8T^ V the free 
perturbation 

From the physical point of view, the intrinsic perturba- 
tion ST^ U corresponds to the irregularities in the matter an d 
that creates the background, i.e. to deviations of homo- 
geneity and isotropy in the matter which produces the 



5T {i)v = 5T (h)u + ST y,)v 

fJi fj, fj, 



(22) 



where ST™" and 5T^ )v are the pieces linear in the met- 
ric perturbation h^, and in the scalar field perturbation 
Scj) respectively. The covariant form of these pieces is 



8Tl h > = (p + p) 



1 _ 

2 UaU/3 



h 01 ? 8" - u^u a h av 



(23) 



and 



8" 

A* 



dV(<t>) 
d4> 



(24) 



Let us discuss now how the perturbations behave under 
gauge transformations induced by infinitesimal coordi- 
nate transformations. Given an infinitesimal coordinate 
transformation x^ — ► x'^ = x M — £ M (x) , the correspond- 
ing gauge transformations of the metric and scalar field 
perturbations are 



is any solution of the homoge- an( j 



(25) 



(26) 



Then, the Einstein tensor perturbation and the intrinsic 
energy-momentum perturbation transform accordingly 



SG' 



8G u + C$ G " 



(27) 



STW* = STjp" + C K f/ , (28) 
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where is the Lie derivative with respect to the in- 
finitesimal vector field £ M . We have assumed, as usual, 
that the zeroth order parts of all tensors remain invariant 
in a fixed coordinate system on the background space- 
time manifold (for example (rj,x), but not necessarily 
this one), and that the effect of infinitesimal coordinate 
transformations is fully charged to the perturbations of 
tensors. Also notice that the transformation laws j27| ) 
and ( p8| ) , can be directly checked by replacing the trans- 
formation laws ( p5f ) and (26) into the functional expres- 
sions (|l|), @ and @. 

It is important to observe that the free perturbation 
8T^ V is gauge invariant. This is compulsory because 
the total perturbation ST" has the same transforma- 
tion law j28| ) as the intrinsic perturbation 8T^ . This 
happens due to the fact that the gauge transformation 
of 8T^ V is of second order in perturbations, since this 
energy-momentum tensor has no zeroth order counter- 
part. Moreover, the zeroth order Einstein equations ( |l6| ) 
imply that the combination 



8G^ = SG^ 



8tt£ ST^ 



(29) 



is gauge invariant. Therefore, by sending the intrinsic 
perturbation to the l.h.s. of (jl7|), we obtain the equations 



(30) 



in which both sides are gauge invariant. Thus equations 
( |30| ) for cosmological perturbations are covariant with 
respect to finite general coordinate transformations, and 
gauge invariant with respect to infinitesimal coordinate 
changes around a fixed coordinate system on the back- 
ground space-time. 

Also notice that both sides of (|3(]) , are covariantly con- 
served with respect to the background metric. 



the gauge. Although all gauges contain the same physics, 
most of them contain it in a complicated way, meaning 
that the independent functions which solve the equations 
that remain after the gauge fixing, have spurious compli- 
cations related to the gauge, but not related to real com- 
plications in the physics. The situation here is similar 
to that of the simple system of gauge invariant Maxwell 
equations for the photon field OA^ — d^d^A") = J^. 
In Lorentz gauge [d v A v = 0), one has nice solutions for 
these equations, but one could complicate the functional 
form of the solutions by choosing an inappropriate (but 
possible) gauge fixing condition. 

In the case of the equations for cosmological perturba- 
tions (^), the appropriate gauge simplifying the equa- 
tions is not so obvious as in the case of electrodynamics. 
The strategy that we follow is inspired by the way in 
which harmonic gauge works to simplify the equations 
for perturbations around Minkowski space-time. First, 
we set the covariant gauge fixing condition 



fl \v 



B„ 



(31) 



where the field is yet unspecified. Then, if does 
not depend on the derivatives of the metric perturba- 
tion ij)^ |„, this condition eliminates all terms in second 
derivatives coming from the second, third and fourth 
terms in the expression (|l^) for 8G^ . Thus, only the 
flat dAlembertian remains as second order differential 



operator in the left hand side of (30). In addition, work- 
ing in coordinates {f],x), we fix the B^ field to simplify 
the equations ( J30| ) as much as possible. After the long 
but straightforward exercise of writing 8G ^ given by the 
covariant expressions ([l8]), ((2^), ( |24| ) and (|2^) in coordi- 
nates (77, x), the analysis of the resulting expression shows 
that a great simplification is achieved by the choice 



4. Choosing the gauge and decoupling of the 
equations 

Now we are finally going to give the covariant and 
gauge invariant equations ( |30| ) a nice definite form. This 
entails a combination of two things: choosing a conve- 
nient coordinate system on the background, and fixing 



B. =-2Hu v i) u 



(32) 



which has also been written in covariant form. Then, 
using the choice (|32|) to specify the gauge ([u]), the cos- 
mological perturbations equations ( p0|) in this gauge, and 
in coordinates (rj,x), take the simple form 



J 



□ V/ - 2 d v ^» + 2n [81 8" ^-81 V - 8 V Q W°] + 32ttS 5° 8% 



= -WttG a 2 8Tjf> 



(33) 



r 



where □ = 
dAlembertian. 



-d 2 /drj 2 + V 2 is the Minkowski 



In addition the covariant conservation of 8T^ V , when 
applied to equations (B3|) yields an equation for the per- 



turbation 8(j), which can also be obtained by perturbing 
the full equation of motion for the field <\>. This last equa- 
tion can be skipped since it is implicit in the equations 
(|33|), as can be checked by direct computation and tak- 
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ing into account the gauge fixing condition. Instead the 
scalar field perturbation 5<p can be expressed in terms of 
the metric perturbations through the gauge fixing condi- 
tion. In coordinates (r],x), we obtain 



be recast in the form 



□ + n 2 + o 



(oV/) = -WnGaHT^ , (35) 



8<j) 



1 



IGttQ i 



^ "+ n(2Vg- V) (34) 



□ + ii 2 - q 



Finally, replacing (M) in (p3|) , these ten equations can and 



aip L 



16tt£ a 3 5T { F)J , (36) 



□ 





^(2Vo°-^) 



16tt0 — (2 £T (f )0 - <ST (F) ) - {d^/ + d^j 



(37) 



Equations (p5[), ( p6| ) and fl37j ) are a decoupled system of 
hyperbolic partial differential equations for the ten com- 
ponents of the metric perturbation. They can be exactly 
integrated and their solutions expressed as retarded po- 
tentials using the retarded Green's functions for the three 
differential operators appearing on the left hand sides. 
Equations ( p5| ) and (|36|), for the ipl and V'o components, 
must be solved first because some combinations of the 
spatial components of the metric perturbation enter as a 
source in equation (p7|). 



The decoupled system of hyperbolic equations (35), 
( |36| ) and (^) is the central result of this paper. They 
give a unified description of all three types of metric per- 
turbations: scalar, vector and tensor, and they have a 
high degree of generality since they are valid for any spa- 
tially flat FRW background and any perturbation that 
can be decomposed as a scalar field perturbation plus an 
arbitrary (weak) perturbation corresponding to geodesi- 
cally moving matter on the background. The key ingre- 
dient in obtaining this system has been the appropriate 
selection of the gauge fixing condition given by ( |3l| ) and 
( |32| ) . The residual gauge invariance remaining under this 
gauge fixing condition has been imposed is given by the 
vector fields £ p fulfilling the equation 



In the case in which the free perturbation ST^ V van- 
ishes, the gauge fixing condition plus the residual gauge 
invariance reduce the number of physical degrees of free- 
dom contained in the metric perturbation tj) u plus the 
scalar field perturbation 8<j>, down two three: a scalar per- 
turbation plus two polarizations for gravitational waves 
on the background. 



5. Universes with a linear equation of state 

To obtain a definite expression for the retarded Green's 
functions solving equations (35), fl36| ) and (|37j), a partic- 
ular background model has to be chosen. In the case of 
universes obeying a linear equation of state p = ap, these 
equations take a particular simple form, and in addition 
the operators on the left hand side of equations ( |36| ) and 
([57]) coincide. Thus, in this case only two Green's func- 
tions for the two differential operators 



CI 2 



□ 



2 -6a 1 
(1 + 3a) 2 rf 



(41) 



and 



□£„ + 2 n d n ^ + 2n - 5° £o) = , (38) 
that can be decomposed in the two simple equations 



and 



□ 



h 2 - n 



si 2 + n 



K) = 



(of) = , 



(39) 



(40) 



involving the same differential operators (and hence the 
same mode solutions) as d3fj) and 



n 2 



n 



□ 



6 + 6a 1 
(1 + 3a) 2 rf 



(42) 



are required. 

The Green's functions for these operators can be ex- 
pressed as superposition of homogeneous mode solutions, 
which are given by Bessel functions of indices v± — 
±(3-3a)/(2+6a)and^ 2 = ±(3a+5)/(2+6a). Moreover 
in the particular cases of de Sitter (a = —1) [pl|, dust 
(a = 0) |6), and radiation (a = 1/3) dominated cosmo- 
logical models, the Bessel indices are the simplest half in- 
tegers (u! = ±3/2,^ 2 = ±1/2), ( Vl = ±3/2,i/ 2 = ±5/2), 
and [y\ = ±l/2,i/2 = ±3/2) respectively. 
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Thus, to integrate the perturbations equations for de 
Sitter, dust, and radiation backgrounds we only require 
the three retarded Green's functions solving 



□ + G^\x,x') = -S^(x-x') 



(43) 



for A = 0, 2, 6 and where x = (rj, x). 

Then G^(x,x') is the well known retarded Green's 

function for Minkowski, and G^\x,x') and (x, x') 
can be obtained by elementary QFT methods as ex- 
plained in |15|] , with the result 



Gp,\x, x') 



1 



4-7T \ x — X' 



' Hv ~ V ~ \x — x'\ 



1 



+ a lO{r,- V ' -\x-x'\) , (44) 

47T 77/7' 



and 



G R (x,x') — 



1 



47T \x - 

3 



1 



' ^iv — i]' — \x — x'\) 
1 \x~x'\ 2 



Thus, G^\x,x') and (x, x') are given by the 
Minkowski retarded Green's function with support on the 
past light cone plus an additional term with support in 
the interior of the past light cone. Then the metric per- 
turbation tp v can be expressed as retarded potentials by 
means of these Green's functions. 



7f V 

x9{j] — r\ — \x — x'\) 
(6), 



rfrj' 



(45) 



6. Concluding renarks 

We have analyzed cosmological perturbations on a spa- 
tially flat FRW background given by a scalar field per- 
turbation plus an arbitrary (weak) perturbation corre- 
sponding to test matter moving geodesically on this back- 
ground. We have shown that by an apropriate choice of 
the gauge the cosmological perturbations obey decoupled 
hyperbolic equations. The class of perturbations consid- 
ered have a high degree of generality, although it is not 
fully general since a covariantly conserved tensor, lin- 
early dependent in the metric perturbation, could still 
be added to the perturbing energy-momentum tensor. 

On the other hand since the treatment and the gauge 
fixing condition that we have introduced is covariant with 
respect to finite coordinate transformations on the back- 
ground, it provides a very suitable framework to develop 
linear quantum gravity around a spatially flat FRW back- 
ground. 

It is important to fully elucidate the relationships of 
the present treatment with the standard Lifshitz-Bardeen 
treatment in terms of gauge invariant perturbations and 
the decomposition of the metric perturbation in scalar, 
vector and tensor parts. This task together with more 
details of calulations and applications are being worked 
out and will be the subject of a forthcomming paper. 
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